Interaction induced dimerization in zigzag single wall carbon nanotubes 
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We derive a low-energy effective model of metallic zigzag carbon nanotubes at half filling. We show 
that there are three important features characterizing the low-energy properties of these systems: 
the long-range Coulomb interaction, umklapp scattering and an explicit dimerization generated by 
interactions. The ratio of the dimerization induced gap and the Mott gap induced by the umklapp 
interactions is dependent on the radius of the nanotube and can drive the system through a quantum 
phase transition with SU(2)i quantum symmetry. We consider the physical properties of the phases 
on either side of this transition which should be relevant for realistic nanotubes. 

PACS numbers: 71.10.Pm, 73.63.Fg 



I. INTRODUCTION 

Since their discovery? 1 - carbon nanotubes have at- 
tracted a great amount of attention, both theoretically 
and experimentally J 2 - It is widely appreciated that single- 
wall nanotubes (SWNT) constitute almost ideal systems 
where peculiar effects specific to strong correlations in 
one dimension (ID) can be observed. They are struc- 
turally ID objects built by wrapping a sheet of graphene 
into a cylinder. The type of wrapping is characterized 
by the superlattice vector (n, to). From the point of view 
of strong correlations metallic nanotubes are of particu- 
lar interest. In the absence of interactions between the 
electrons? 3 - the band structure is indeed metallic for all 
{n, n) (or armchair) nanotubes, as well as for the (n, —n) 
(or zigzag) nanotubes with n being multiple of 3. In all 
other cases there is band gap which at large n is of the 
order of (10/n)eV. 

The low-energy effective theory for correlated metal- 
lic SWNTs away from half-filling 4 shows that the long- 
range Coulomb interaction converts the nanotubes into 
Tomonaga-Luttinger liquids implying various scaling 
laws for conductance, which have since been verified 
experimentally^ At half-filling the effects of strong cor- 
relations become even more pronounced in the presence 
of the poorly screened Coulomb interaction. The lat- 
ter makes the umklapp scattering processes strongly rel- 
evant giving rise to sizable Mott gaps£ The underlying 
strong-coupling phases for armchair nanotubes were re- 
cently classified in Ref. 0. 

In this paper, we study the effects of the long-range 
Coulomb interaction in half-filled zigzag SWNTs (n divis- 
ible by 3). The new feature that makes this case qualita- 
tively different from the armchair nanotubes is an explicit 
dimerization that originates from the nearest-neighbor 
interaction and gives rise to a single-particle gap of the 
order of 1/n. Notice that curvature effects^ which are 
known to lead to much smaller band gaps, ~ l/^ 2 , can 
therefore be safely neglected. We will explain how the 
dimerization comes about, discuss its interplay with the 



Mott gaps and show that this competition can result in a 
quantum phase transition similar to the one studied pre- 
viously in the context of dimerized spin ladders j^iiil Mak- 
ing reasonable approximations for the relative strengths 
of different interaction terms for a zigzag SWNT shows 
that the system can occur very close to the quantum crit- 
icality, and may even be tuned to reach it by changing 
the radius of the nanotube; however a realistic estima- 
tion of the parameters by either experiment or numerical 
techniques is beyond the scope of this paper. 

Below we will follow the same strategy as in the arm- 
chair casei paying particular attention to the differences. 
The structure of the paper is as follows: First we consider 
the band-structure of the zigzag SWNT within a tight- 
binding model and show that, in the low-energy limit, 
the latter is equivalent to an effective two-chain model. 
We then very carefully consider the structure of the in- 
teraction, taking into account both the long-range tail 
of the unscreened Coulomb interaction and the details 
of the short-range component of the interaction which 
depend in an important way on the lattice. We then 
solve this model by using bosonization, employing the 
adiabatic approximation to treat neutral collective exci- 
tations, and then concentrating on energies well below 
the charge gap, by refermionizing the remaining theory. 
This allows us to identify the nature of the phase tran- 
sition generated by the dimerization term, and extract 
some physical properties of such nanotubes. 



II. MAPPING ONTO LOW ENERGY 
EFFECTIVE THEORY 

In this section, we single out those two bands that cross 
the Fermi level in an undoped isolated zigzag SWNT, 
which allows us to describe the relevant part of the 
spectrum in terms of an equivalent two-chain fermionic 
model. We carefully analyze the details of the lattice 
structure that determine the peculiarities of the interac- 
tion terms. We then pass to the continuum limit and em- 
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FIG. 1: Real space structure of a zigzag nanotube. In this 
paper the chain direction is chosen to be x, meaning that 
the axis of the graphite sheet have been rotated by 90° as 
compared to the references that deal with the armchair case. 



ploy the bosonization technique to derive the low-energy 
effective field theory for a metallic zigzag SWNT. 



Kinetic term and two chain model in zigzag 
nanotubes 



We begin with a tight-binding model on the honey- 
comb graphite lattice shown in Fig.l. The structure con- 
sists of an underlying triangular lattice with two inequiv- 
alent Carbon positions, labeled by a and b. It proves 
convenient to double the unit cell such that the new cell 
contains four atoms, a, b, a, b. In this picture the unit cell 
is rectangular, so the two momenta, k x and k y , can be 
considered as independent. Wrapping a (N x n) graphite 
sheet into a (n, —n) zigzag nanotube quantizes the elec- 
tron momentum in the direction of the superlattice vector 
C = n(ai — a2), i.e. the y direction in our Cartesian axis. 
Carrying out a partial Fourier transform in this direction 
gives a Hamiltonian describing n ID bands (see e.g. Ref. 
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5o = -EEM') [cl(l)c q (l + l) + h.c] 



(1) 



Here I labels the sites of an effective ID lattice along the 
tube's axis, the lattice spacing being xi+\ — xi = b = 
\/3/4, q is the band index taking n integer values (\q\ < 
ri/2), and t q (l) = to(q) + (—l) l A(q) is an alternating 
hopping amplitude whose uniform and staggered parts 
are given by 



to(q) 
Ml) 



<H + 2t± cos(irq/n)] 
t|i — 2t± cos(nq/n)] 



(2) 
(3) 



Thus the spectrum is generically gapped, the gap of 
the qth band being determined by the dimerization am- 



plitude ([3]). However, for two bands labeled by q = ±Q = 
±n/3 (n = 3m) the gap is minimal, proportional to the 
difference t\\ — t±, and is solely due to curvature effects. 
This difference was shown£ to be of the order 1/n 2 and 
is much smaller than the gaps of all the other bands. For 
this reason, and also in view of the fact that the interac- 
tions are of the order 1/n or larger, one can set t\\ = ij_ 
and thus make the lowest-energy part of the spectrum 
described in terms of two decoupled, translationally in- 
variant chains with a simple half- filled cosine band, via 
the mapping 



b{m x ,m y ) 
a(m x ,m y ) 
b(m x ,m y ) 



a(m x ,m y ) — > -= 



l j= [e iQv c + (Am x )+e- iQy C-{Am x )] 
In 

\= [e iQv c+(4m x + 1) + e- iQv c-{4m x + 1)] 
e iQ(!/+i/2) c+ (4 mi + 2) 

+e- iQto+1/2) c_(4m x + 2) 

e iQ(y+1/2) c + {4m x + 3) 
+fr iQto+i/2) c _ (4ma; + 3) _ (4) 



This picture is in contrast with the case of armchair nan- 
otubes where the kinetic part of the low-energy theory 
is that of two strongly coupled chains (ladder). Even 
though curvature effects can be neglected, the dimeriza- 
tion A should be retained in the Hamiltonian. This fol- 
lows from the fact that interactions generate an identi- 
cal term in the low-energy effective theory, but with a 
much larger amplitude, ~ 1/n. The next two sections 
will demonstrate this phenomena, which is one of the 
central results of this paper. 



B. Interactions 

There are two important contributions to the inter- 
actions. Firstly, in an isolated nanotube, the Coulomb 
interaction is unscreened, so one must take into account 
the effects of its long-range tail. Secondly, at half-filling, 
Umklapp processes originating from the short-range part 
of the Coulomb interaction play a crucial role. Below we 
consider these two contributions separately. 

The long-range part of the interaction is insensitive to 
the details of the lattice and is simply given by 



H Co ui = Y,n(l)U{l-l')n(l'), 



(5) 



where n(l) = V CT 



is the total density on lattice 



site I in the two chain mapping, and U{x) ~ e 2 /\x\. 

The effects of the short-range part of the interaction 
are best elucidated by considering the minimal model 
that captures all essential physics of the problem. It turns 
out that a two-parameter model with on-site (U) and 
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FIG. 2: A cartoon representation of the origin of the staggered 
interaction 



nearest-neighbour (V) couplings fully represents the most 
general structure of local interaction in the low-energy 
limit, 

Hi = U 22 n^(m x ,m v )ni(m x ,my) 

m x ,m y 

+ V X! n { m x,m y )n(m' xl m' y ). (6) 



Here m x ,m y label the carbon atoms on the original 2D 
lattice, the summation in the second term goes over 
nearest neighbours, and n(m x ,m y ) — n^(m x ,m y ) + 
ni(n x , n y ) is the electron density on site (m x , m y ). When 
the V-term is projected onto the two low-energy chains, 
we see an interesting effect as illustrated in Fig. [2j each 
atom a (a) interacts with one atom b (b) and with two 
atoms b (b). In other words, when mapping onto the 
two-chain low-energy sector, the interaction V a b = V-^ ^ 
Vt = Vab- This means that in the low-energy theory 
the V interaction will contain not only a uniform part 
but also a staggered part. The latter is entirely due to 
the way we break the C3 symmetry of the 2D lattice by 
wrapping it to form a zigzag nanotube; the effect is not 
present in the armchair case. Carrying out the mapping 
(|3|) to the effective two-chain model, with index fi = ± 
labeling the two chains, yields three terms, 



I ft 



W 
2^ 



H * = -iD-^Wntf + l) 



2n 
V 



(J) 



Ifiaa' 



where Hu is the contribution from the on-site interac- 
tion, Hy is the smooth contribution from the nearest 
neighbor interaction, and H s is the staggered part of the 
nearest neighbor interaction. We now proceed to study 
the ground state and elementary excitations of such an 
interacting two chain model by employing the bosoniza- 
tion technique'. 



C. Chiral Decomposition and Bosonization 

Passing to the continuum limit, we adopt the standard 
description in terms of chiral (right/left) fermions. At 
half- filling kp = ±ir/2b, so that 



£V(0 -> Vb [i l R^{x) + {-ifL^x)] , 



(8) 



where a =|,| is the spin index. This yields H = 
J dx 7io(x), where 



Ho = -iv F (R^dxR^c, - L^ a d x L p 
Cdim = i {w^Lfu, ~ L' R^) . 



23A0dii 



(9) 



To treat interactions nonperturbatively, we employ 
Abelian bosonization based on the correspondence 
R{L)na {lira)- 1 / 2 exp [- iy^r ($ M<7 T 6 M(T )] , where 
$^0- and 0^^ are a pair of mutually dual scalar fields, 
and a is a short z distance cutoff of the bosonic theory^ 
Following Ref. [H, we pass to linear combinations of the 
bosonic fields describing the total and relative charge and 
spin excitations, 



i($ +T +$ +i ±$_ T ±$_i), 



The kinetic energy of chiral Fermions then becomes the 



sum of four Gaussian models (a 



H kin - (VP/2) [(d x Qaf + (d x <S>aY 



(10) 



while the dimerization operator takes the form: 
Odim = (4/7ra) 

x [cos(\/7r$c ) cos(V7T$c ) cos( ) cos(v / 7r3>j") 
+ sin(v / 7T*^ ) sin( % /7r$ f r) sin( % /7T$~) sin(v / 7T$j")]. 

(11) 

As in the armchair case, the long-range part of the 
unscreened Coulomb interaction only involves the total 
charge field 



H, 



Coui = / dx / dy : : . 

t J J \x-y\ 



(12) 
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FIG. 3: The umklapp interaction g-ology: (a) in-chain, (b) 
inter-chain, no spin exchange, (c) inter-chain, spin exchange. 
The solid and dashed lines describe particles with opposite 
chiralities. 



and greatly reduces the scaling dimension of the operator 
cos(v / 47r<i>+). At half filling, Umklapp processes contain 
this operator, and therefore become strongly relevant giv- 
ing rise to a significant increase of all the gaps induced 
by interaction.— i 

Turning to the uniform part of the short-range in- 
teraction ([7]), we single out contributions of umklapp 
processes, all of them proportional to the operator 
cos(v^f7r < 5 ) ( t). This selection can be done by directly 
bosonizing [[7]). and can also be understood physically 
by using the g-ology approach (see Fig. [3]). The result 
is: 

Humki = —r- ^-cos(\/47$+) \g c cos(V^<5>~) 

+ (53 - 5i) cos(\/47r$+) +gi cos(V4tt$7) 
- 53 cos(\/4^e;)}. (13) 

Within our U-V model, the parameters are gi — g c — 
(U — 3V)/n, 53 = U/n. The slightly strange looking 
structure in the spin sector is an artifact of treating 
the non-Abelian SU(2) symmetry group by means of 
Abelian bosonization; this symmetry becomes manifest 
upon refermionization (see next section). 

Appendix [X] contains a discussion of the full structure 
of the interaction including the less relevant terms which 
we will ignore in the remainder of this paper. The only 
feature generated by this extra complexity is a minor log- 
arithmic renormalization of the masses generated in the 
theory. This may be important for estimating quantities 
to be measured experimentally, but will not make a dif- 
ference to any of the universal properties of the gapped 
ground state phases of the system and the quantum crit- 
icality separating them. 

The staggered part of interaction Q is contributed 
both by the spin-exchange and non-spin-exchange scat- 
tering processes. In the effective one-dimensional picture, 
the non-spin-exchange process is proportional to the sum 
l) l n{l)n{l + 1). Here n(l) is the total local den- 
sity that transforms in the continuum limit to b[J(x) + 
(-l) l M(x)}, where J = £ M(T (■ R^R^ ■ + : L^L^ :) 
and M — ^^{R^L^a + h.c). The above sum then be- 
comes b J dx [M, J] x ,x+a- The emerging point-split com- 
mutator is then estimated using Operator Product Ex- 
pansion (see Appendix [B]) and is proportional to the 
single-particle dimerization operator 

[M, J] x/x+a = lim {M (x)J(x + a) - J(x)M (x + a)} 



= 20 dim {x). (14) 

A similar procedure results in the same operator for the 
spin-exchange processes. Within the U-V model the 
dimerization amplitude is gA = V/2n. We must stress 
again here that although Odim is a single particle dimer- 
ization operator, g& is an effective coupling constant, 
arising from projecting the interactions on the hexagonal 
lattice onto the low-energy theory relevant for a zigzag 
nanotube. The bare coupling constant associated with 
this operator arising from curvature effects may then be 
safely neglected, as it is of order 1/n 2 , much less than 
the effective gA- 

Thus we arrive at the effective model describing the 
zigzag nanotube at half filling in the scaling limit, 

H = Wkin + HqouI + H um kl + , (15) 

where Ha = gA Chim- 
in. ADIABATIC APPROXIMATION AND 
MAPPING TO SPIN SECTOR 

The unscreened Coulomb interaction strongly en- 
hances the velocity of the symmetric charge mode rel- 
ative to the group velocities of collective excitations in 
other channels*^ Hence, from the point of view of the 
field <&+ , all other degrees of freedom can be regarded as 
static. This allows one to employ an adiabatic approxi- 
mation and obtain the low-energy dynamics of the model 
by integrating out the total charge mode. Such a proce- 
dure has already been discussed for narrow-gap SWNTs 
away from half-filling 15 (7i U mki = 0) and for a half-filled 
armchair nanotube 7 (Ha = 0). The general analysis of 
the model (fT5*)) represents a more complicated task due 
to the presence of two strongly relevant perturbations 
- the Umklapp and dimerization terms, each containing 
the field $+. 

In this paper, we will specialize to the limit where 
W um ki is the principal term, responsible for the forma- 
tion of a fully gapped Mott phase, and then address the 
effect of Ha as a perturbation on this phase. 

Such an approach can be justified in our U-V model 
by choosing U ^> V > 0, because the structure of the 
Umklapp interaction (|13p contains parts proportional to 
U, whereas the dimerization includes only the smaller V. 
The details of integrating out $+ are well explained in 
Ref. 0- Here we give a basic physical argument that 
leads us to the correct result. By examining Eq. (fl~3f and 
making the approximation that all other fields are static 
during the typical fluctuation time of <f>J, we obtain a 
modified sine-Gordon model describing the total charge 
field, 

H ~ fto[$+] - AcosV4^$+. (16) 

Here Ho[&~^} is the sum of a standard Gaussian model 
and the nonlocal Coulomb term (fl"2|) . whereas A is a com- 
bination of the slower fields which, within the adiabatic 
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approximation, can be replaced by a constant. Due to 
the high velocity and low effective scaling dimension of 
the field <f>+, all the spectral gaps (i.e. those of solitons 
and breathers) in the total charge sector are very large. 
Therefore, to examine the low-energy dynamics of neu- 
tral modes, one can consider the field to be locked at 
<3>+ = and replace cosines of this field by their expec- 
tation values. 

It is important at this stage to consider the physics we 
are discarding by making this approximation. Charged 
solitons in the field do indeed exist in the theory. 
Furthermore, electron like quasi-particles described by a 
half soliton simultaneously in each of the sectors of the 
theory are also present. Finally, there will also be plas- 
monic excitations corresponding to breathers modes in 
Hamiltonian (|16p . However, due to the strong renormal- 
ization of the field, each of these excitations will have 
a large mass. 

Thus, as we are interested in the low-energy properties 
of the model, the only excitations that remain are neutral 
collective modes: the relative charge "vortex"— modes, 
and the spin degrees of freedom (which as we will see 
later can be separated into a singlet mode and a triplet 
mode), and it is in this region of energies below any of 
the charged excitations that the adiabatic approximation 
(and therefore the description in terms only of neutral 
modes) is valid. This also explains why the dimcrization 
term Ha with a magnitude much less than any of the 
charged excitations can not turn the system into a trivial 
band-insulator: everything occurs within the Mott phase. 
We now proceed to analyse the low-energy phase diagram 
of such a model, and show the existence of an interesting 
non-trivial quantum phase transition. 



A. The Umklapp term 

Following the preocedure of Ref. 0, one can identify the 
remaining neutral collective modes at g& = by intro- 
ducing four real (Majorana) fermions \r l (* = 0' 1> 2> 3) 
and refermionizing the theory. The resulting Hamilto- 
nian displays an SU(2) x Z2 symmetry in the spin sector: 



m f 



COS V 47T$ 



- ^("Xk^xk + Xl^xXl) " im t XflX2 

a=l 

(17) 



i=0 



imsXRXL 



Here mf = (cos(V47r$+))(7i/7ra is the mass of the 
relative charge ("vortex") excitation, whereas mt = 
(cos(V47r$+))( i gi- i g 3 )/7rQ! and m s = (cos(\/47r<I>+))(#3 + 
gi)/na are the masses of the spin-triplet and spin-singlet 
modes, respectively. In terms of our U — V model, the 
masses are parametrized as follows: 



C 



2U-3V 



-C 



where C = (cos(\/47r$+)) w 1. 

Associated with the Majorana fermions are four Ising 
models characterized by order and disorder parameters, 
<Ji and /li. We use the correspondence: 



cos Vn&f = 
cosy^O^ = \i\02 
cos \pK§~ — a (T3 
cos v^Q^ = Mo 03 



sin Vn@t = 

sin = MoM3 

sinv^i'S.r = coM3- 



(19) 



The sign of the Majorana mass indicates whether the 
corresponding Ising model is ordered (m < 0, (a) =/= 0, 
(fi) = 0) or disordered (m > 0, ((i) ^ 0, (a) = 0)±&. 

The fermionic part of the Hamiltonian (|1T[) has the 
same structure as that of a two-leg antifcrromagnetic 
spin-1/2 ladder^ Indeed, assuming that U/V ^> 1, we 
find that m t < and m s > so that (a a ) 7^ for a 
1.2. .'i and (/io) 7^ 0. Furthermore, rrif > meaning that 
(cos(v / 7r ( E ) t r)) 7^ 0. These signs of the masses indicate a 
spin-liquid behavior of the systemi 16 ' 17 We can further 
quantify this by defining the staggered magnetization as 
a suitably averaged difference between the local spin den- 
sities on the two sublattices, n (r) = S (r) — S(,(r). Pro- 
jecting this onto the low-energy sector of the model gives 

n ~ cos(^/ir$ c - )fjt {^i<7 2 cr 3 , o^o^ eric^Ma)- (20) 

The two-point correlation function of n - displays a co- 
herent magnon peak with a mass gap \mt\- 



B. Effect of dimerization term 

After refermionization the dimerization operator ac- 
quires the following low-energy form: 



3( V / 7T$ C )AiOMlM2A«3- 



(21) 



In the spin- liquid phase formed by the umklapp processes 
this operator can be further simplified by replacing in the 
leading order the operators cos(- s /¥<i> f r) and /j,q by their 
expectation values. Thus, projecting 0<jim onto the spin- 
triplet sector yields 

= h e Q ^i/X2M3, 
4<?A 



^eff = 



(ira) 2 



cos(v^$ c + ))(cos(V^<f c "))(^o). (22) 



(18) 



We see that the dimerization term competes with umk- 
lapp processes which support the ground state with 
{a t ) ± 0, = 0. 

The resulting theory thus involves only the spin-triplet 
degrees of freedom and actually coincides with the prob- 
lem of a two-leg spin ladder with explicit dimerization, 9 
or, equivalently, a weakly dimerized spin-1 chain with 
a small Haldane gapJ^ It was shown that increasing 
the dimerization can drive the system towards a quan- 
tum criticality belonging to the universality class of the 
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FIG. 4: Average dimerization as a function of g = m&/mt 



at criticality 
slight deviation 
from criticality 



larity in the derivative at the critical point g = g c . 



being nonzero in both gapped phases, displays a singu- 



FIG. 5: C v /T as a function of T. 
the behavior exactly at criticality g 

slightly away from criticality. The crossover temperature is 
given by the smallest gap in the system m ~ \g — g, 



The solid line denotes 
= g c , the dashed line is 



|2/3 



SU(2)i Wess-Zumino-Novikov-Witten model with cen- 
tral charge 1. Close to criticality the universal lowest- 
energy properties of the system are those of a single, ex- 
plicitly dimerized, antiferromagnetic spin- 1/2 chain, with 
the dimerization changing its sign across the transition. 

The location of the critical point can be estimated by 
requiring that the mass gap generated by the dimeriza- 
tion term alone and the triplet mass become of the same 
order. Using standard scaling arguments we find that in 
Eq. (|22[) h e Q ~ |m/| 1 / 4 |m s | 1/ ' 8 s , A, where all parameters 
in the right-hand side are proportional to 1/n. Therefore 
the dimerization gap scales as 

m dim ~ h*J£ a oc (1/n) 11 / 13 , (23) 

and so, within our U V approximation 




where U — Ua/v and V — Ua/v are dimensionless cou- 
pling constants. 

This estimate shows that for sufficiently large radius 
n the dimerization is always dominant and the system 
occurs in a gapped phase whose properties are governed 
by the operator Odim- Whether a crossover to the spin- 
liquid phase through the SU(2)i criticality can occur as 
a function of the nanotube's radius depends on various 
nonuniversal prefactors. As shown in Appendix an 
improved estimate of the critical line can be obtained by 
taking less relevant (non-umklapp) terms into account. 
The ratio of the dimerization and triplet mass gaps then 
becomes 

m dim [V/Uf' 13 
\m t \ (U/n) 15 / 13 \n{n/U)' 1 ' 

The new estimate shows that the quantum transition sce- 
nario taking place upon decreasing n at large U is a real 
possibility. 



C. Physical Observables 

The ratio g — m^/vnt acts as a control parameter 
for the quantum phase transition, with a critical point 
g c ~ 1. As we vary g, the dimerization itself shows a 
singularity in its derivative as one passes through the 
transition point itself^ according to 

(Odim) = (Odim(g = g c ))+a\g-g c \ 1/3 sign(g-g c ), (26) 

where a is some constant. This is qualitatively shown in 

Fig. a 

The parameter g can be varied by varying the radius 
of the nanotube, n, or by changing the interaction cou- 
pling constants; the latter can be achieved by, for exam- 
ple, applying pressure or stretching the nanotube. The 
dimerization of the nanotube will show up as a Peierls dis- 
tortion of the lattice, which should be measurable with 
low-temperature STM. 

The proximity to the quantum critical point will also 
show up in thermodynamic quantities, such as specific 
heat. Exactly at the critical point, the specific heat is 
linear in temperature^ 

C v (g = g c ) = c—, (27) 

where c = I is the central charge of the criticality, and v 
is the Fermi velocity. Away from criticality, the specific 
heat is exponentially small at sufficiently low tempera- 
ture T < m, 

G v ~ T- 3 ' 2 e~ m / T , (28) 

where m is the smallest gap in the system, which near 
criticality is giveniS by m ~ \g — g c \ 2 ^ 3 - The specific heat 
crosses over to a linear dependence on temperature when 
T > m. This is plotted in Figure [5] 

IV. CONCLUSIONS 

We have shown that in zigzag carbon nanotubes 
electron-electron interactions play an important role in 
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the ground state of the system. Firstly, there is the long- 
range unscreened Coulomb interaction, which sets a hier- 
archy of energy scales in the problem. Secondly, there are 
umklapp processes (in the undoped case) which gap all 
collective excitations in the nanotube. Finally, there is 
an explicit dimerization in the interaction, scaling as 1/n 
due to the way that we break the C3 symmetry when we 
wrap the zigzag nanotube. This dimerization originates 
from interactions and affects the effective low-energy ac- 
tion; there it directly competes with Umklapp processes 
supporting a Mott-like insulating phase with a spin-liquid 
structure of collective excitations. 

As the relative strength of the dimerization is in- 
creased, the system can exhibit a nontrivial quantum 
criticality with an SU (2)i symmetry. On one side of 
this critical point, the system is in a Haldane spin liquid 
phase, on the other side of the critical point, the ground 
state is dimerized. The position in the phase diagram de- 
pends on the radius of the nanotube, with reasonable as- 
sumptions about parameters showing we are always near 
the QCP. Therefore, small external perturbations may be 
able to drive the real system of a zigzag carbon nanotube 
to an exciting theoretical quantum criticality. 
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m f w C— (A2) 
n 

where C — (cos >/4n$>+) rj 1, so we can assume that this 
field <i>~ becomes locked, and the expectation value 

(cos V4n$c ) w to f In — . (A3) 

TO/ 

After refermionization, AH will contribute to both the 
singlet mass and the triplet mass. As the singlet gap is 
already of order U/n, a change of order (U/n) 2 is in- 
significant. However, the triplet mass will become 




We considered in the main text the case where V/n 
(U/n) 2 ln(n/U). If U is large enough, then the opposite 
limit may hold true, in which case the ratio of the dimer- 
ization and the triplet masses modifies and becomes given 
by formula ([25]) . 
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APPENDIX A: FULL STRUCTURE OF 
INTERACTION 



APPENDIX B: THE OPERATOR PRODUCT 
EXPANSION 

Here, we use the Operator Product Expansion (OPE) 
in the fermionic basis to show that the staggered inter- 
action gives rise to a single particle dimerization term in 
the continuous limit. 

The bare Hamiltonian in the Fermionic basis is 



I dx [Rl a d x R^ - Ll„d x L^} , (Bl) 



The long range Coulomb interaction means that the 
gap in the total charge excitations is the largest en- 
ergy scale in the problem, and consequently the umk- 
lapp terms which involve cos\/47r$+ are the most rele- 
vant terms which gap the remaining sectors of the the- 
ory. However, as the term which gaps the triplet sector is 
proportional to V <C U, less relevant terms may strongly 
renormalize the triplet mass. In this section, we continue 
to assume that U ^S> V, so that we only look at correc- 
tions proportional to U. 

Bosonizing the interaction Hjj ([7]) gives rise to the fol- 
lowing addition to the Umklapp Hamiltonian (|13p 



AH = U dx cos V47r$ r cos V 47r$, + cos V 47rO 



(Al) 



where U is a dimensionless interaction constant U 
Ub/^a) 2 . 



so the Matsubara Green's functions are 



2tt(z — w) 



'"" -'"M'^"" 2tt(z-w) 

(B2) 



where 



z = vpT + ix, z = vpT — ix, 



(B3) 



and similarly for w, are the complex coordinates written 
in a manifestly Lorentz invarient way, with r being the 
imaginary time. 

The first staggered piece in Eq. {7]) takes the form 



(-l) l n(l)n(l + 1), 



(B4) 



8 



where n(l) — ^ n^ a (l) is the total electron density at 
site In the continuum limit, this becomes 



b [M{x] J(x + a) - J(x)M(x + a)} , (B5) 



where 



J(x) = + 

m{x) = J2( R U L ^+ L U R 



(B6) 



jiCT 



Therefore, in the expansion of (|B5|) . we will need the 
following OPE: 



: L\z)L(z) :: lS{w)R{w) : 
{L{z)L\w)) : LHz)R(z) := 



1 



LHz)L(z) :: B*(w)L(w) : 
-(L^(z)L(w)) : RHz)L(z) := 



2ir{z - w) 
-1 



2tt{z — w) 



L^z)R(z) :, 
: R){z)L{z) 



(B7) 



where the flavor and spin subscripts must all be identi- 
cal. Similar OPE's hold for the antianalytic part of the 
expansion, so that 



J(z)M(w) ~ -L / ' 



1 



2ir \z — w z — w / 
x ^ : - L^i?^) : . (B8) 

fX<7 

We recognise the operator as the dimerization operator, 
hence putting r = and evaluating the point-split com- 
mutator in (IB5I). we see that 



[M, J] x .x+a = lim {M(x)J(x + a) - J(x)M(x + a)} 

a— >0 



20dim(x). 



(B9) 



The other component of the staggered interaction 

(-l) l cl a (l)c^ a (l)cl^,(l + l)c^(l + l) 



(BIO) 



can be treated in an identical way. This too turns into the dimerization operator, giving the final result quoted in the 
main text. I 
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